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ABSTRACT 

We  report  new  infrared  radial  velocity  measurements  obtained  with  CSHELL  at  NASA’s  Infrared  Telescope  Facility 
that  reveal  the  M3. 5  dwarf  GJ  867B  to  be  a  single-lined  spectroscopic  binary  with  a  period  of  1.795  ±0.0 17  days.  Its 
velocity  semi-amplitude  of  2 1.4  ±0.5  km  s_1  corresponds  to  a  minimum  mass  of  61  ±  7  Mjup;  the  new  companion, 
which  we  call  GJ  867D,  could  be  a  brown  dwarf.  Stable  astrometric  measurements  of  GJ  867BD  obtained  with 
CTIO’s  0.9  m  telescope  over  the  last  decade  exclude  the  presence  of  any  massive  planetary  companions  (7-18  MJUP) 
with  longer  orbital  periods  (2-10  yr)  for  the  majority  of  orientations.  These  complementary  observations  are  also 
used  to  determine  the  trigonometric  distance  and  proper  motion  of  GJ  867BD;  the  measurements  are  consistent  with 
the  HIPPARCOS  measured  values  of  the  M2  dwarf  GJ  867 AC,  which  is  itself  a  4.1  day  double-lined  spectroscopic 
binary  at  a  projected  separation  of  2  4" 5  (216  AU)  from  GJ  867BD.  These  new  measurements  strengthen  the  case 
that  GJ  867 AC  and  GJ  867BD  are  physically  associated,  making  the  GJ  867  system  one  of  only  four  quadruple 
systems  within  10  pc  of  the  Sun  ( d  =  8.82  ±  0.08  pc)  and  the  only  among  these  with  all  M-dwarf  (or  cooler) 
components. 
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1.  INTRODUCTION 

Although  the  processes  involved  in  single  star  formation  have 
a  wide  variety  of  observational  and  theoretical  constraints  (e.g., 
McKee  &  Ostriker  2007),  the  complex  physics  involved  in  the 
formation  of  binary  and  especially  multiple  star  systems  is  much 
less  understood  (e.g.,  Halbwachs  et  al.  2003;  Tokovinin  2008; 
Duchene  &  Kraus  2013).  Prompt  fragmentation  of  a  collapsing 
prestellar  core  is  favored  as  setting  the  initial  conditions  for 
multiple  star  systems  (e.g.,  Tohline  2002),  though  how  this 
depends  upon  the  initial  conditions  of  the  cloud  core  (e.g., 
metallicity,  mass,  etc.)  is  not  well  understood.  Once  initial 
fragments  form,  subsequent  fragmentation  (possibly  within  a 
disk),  accretion  from  the  remnant  cloud  core,  A'-body  dynamics, 
and  possibly  orbital  migration  have  been  proposed  to  play 
competing  roles  in  setting  the  final  properties  of  multiple  star 
systems  (e.g..  Bate  2000;  Bonnell  2001;  Delgado-Donate  et  al. 
2004;  Sterzik  &  Tokovinin  2002).  While  the  task  of  identifying 
the  relative  importance  of  these  dynamical  phenomenon  remains 
daunting,  useful  constraints  can  nevertheless  be  established  by 
first  determining  fundamental  multiple  star  properties,  such  as 
the  overall  multiplicity  fraction,  the  distribution  of  mass  ratios, 
the  distribution  of  separations,  and  the  mass  dependence  of 
companion  properties.  As  an  example,  one  of  the  most  well- 
established  observational  trends  is  that  the  fraction  of  stars  in 
binary  systems  decreases  toward  lower  stellar  masses;  studies 
of  nearby  stellar  samples  reveal  that  26%  of  M-dwarfs  have 
companions  (Delfosse  et  al.  2004;  also  see  Fischer  &  Marcy 
1992;  Henry  &  McCarthy  1992),  compared  to  46%  for  Sun¬ 
like  stars  (Raghavan  et  al.  2010;  also  see  Duquennoy  &  Mayor 
1991).  A  likely  implication  of  this  is  that  lower  mass  cloud 


cores  fragment  less  often,  despite  any  subsequent  dynamical 
evolution. 

Formation  events  that  lead  to  triple,  quadruple,  and  even 
higher  order  stellar  systems  also  appear  to  occur  less  often 
at  low  stellar  masses.  The  emerging  statistics  are  that  2.1% 
of  M-dwarfs  occur  in  triple  systems,  while  0.03%  occur  in 
quadruple  or  higher  order  systems,  compared  to  5.6%  of 
G-dwarfs  occurring  in  triple  systems  and  2.2%  occurring  in 
quadruple  or  higher  order  systems  (Eggleton  &  Tokovinin 
2008).  While  it  is  tempting  to  interpret  this  as  a  corresponding 
decrease  in  the  frequency  of  multiple  fragmentation  events 
with  lower  cloud  core  mass,  accurate  assessments  of  high-order 
multiplicity  are  much  more  subject  to  observational  biases  and 
incompleteness,  especially  at  low  stellar  masses.  As  an  example 
of  the  incompleteness  that  exists.  Law  et  al.  (2010)  find  a 
surprisingly  high  fraction  of  triple  and  quadruple  systems  among 
a  targeted  sample  of  36  extremely  wide  M-dwarf  pairs,  formerly 
known  only  as  wide  M-dwarf  binary  stars;  the  results  mimic  the 
discovery  that  brown  dwarfs  are  more  often  found  to  be  a  close 
binary  system  when  first  identified  as  a  wide  companion  to  a 
star  (Burgasser  et  al.  2005). 

The  most  direct  way  to  establish  a  more  complete  assessment 
of  the  multiplicity  fraction  over  a  broad  range  of  masses  is  to 
target  volume-limited  samples  of  nearby  stars.  Not  only  are  these 
the  closest,  brightest  stars  with  the  most  precisely  determined 
kinematics,  but  if  trigonometric  distances  are  known,  Malmquist 
biases  can  be  avoided.  Nevertheless,  observational  biases  remain 
even  for  the  closest  stars.  Of  the  348  stars  known  within  10  pc, 
21  are  G-dwarfs  and  239  are  M-dwarfs  (Henry  et  al.  2006). 
But  while  two  of  these  21  G-dwarfs  are  the  primary  stars  in 
quadruple  star  systems  (no  higher  order  G-dwarf  systems  are 
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Table  1 

Quadruple  or  Higher  Order  Systems  within  10  pc 


Primary 

Star 

Spectral  Type 
of  Primary3 

Parallax 

(mas)b 

Configuration0 

GJ  423A 

G0VJ1 

119.51  ±  0.791,  2’  3 

AC-BD 

GJ  570A 

K4V2 

170.62  iO.671-2-  3 

A-BC-D 

GJ  695A 

G5IVJ3 

120.32  ±0.162' 3 

AD-BC 

GJ  867 A 

M2VJ4 

115.01  ±  1.302,  3 

AC-BD 

GJ  644A 

M3.5VJ4 

154.96  ±  0.521,  2-  3-  4 

ABC-D-GJ  643 

Notes. 

“The  letter  “J”  indicates  a  spectral  type  based  on  light  from  a  spatially  unresolved 
pair.  References:  (1 )  Batten  et  al.  1989;  (2)  Gray  et  al.  2006;  (3)  Gliese  &  JahreiB 
1991;  (4)  This  work. 

b  Parallaxes  are  weighted  parallaxes  from  the  following:  (1)  Soderhjelm  (1999); 
(2)  van  Altena  et  al.  (1995);  (3)  van  Leeuwen  (2007);  (4)  Costa  et  al.  (2005); 
(5)  Martin  et  al.  (1998);  (6)  Martin  &  Mignard  (1998). 
c  Widely  separated  pairs  are  indicated  with  a  hyphen. 

known),  only  one  of  the  239  M-dwarfs  is  the  primary  of  a 
quintuple  star  system;  no  M-dwarf  quadruple  star  systems  are 
known  within  10  pc  (see  Table  l).6 

The  deficiency  of  known  higher  order  pairs,  in  part,  stems 
from  a  lack  of  high-dispersion  spectroscopy  that  would  reveal 
short-period  companions.  Early  M-dwarfs  are  faint  in  the 
optical  and  chromospherically  active  making  spectroscopic 
companion  searches  observationally  difficult.  It  becomes  even 
more  challenging  to  search  for  short-period  companions  around 
mid  to  late  M-dwarfs  with  spectroscopy,  as  these  stars  are  even 
fainter  and  tend  to  have  higher  projected  rotational  velocities 
compared  with  earlier  M-dwarfs  (Jenkins  et  al.  2009;  Reiners 
&  Basri  2010;  Reiners  et  al.  2012).  These  biases  have  inhibited 
good  assessments  of  high-order  multiples.  For  example,  of  the 
90  M-dwarfs  included  in  the  spectroscopic  companion  search 
by  Endl  et  al.  (2006),  only  four  stars  have  spectral  types  of  M4 
or  later.  The  HARPS  RV  sample  of  102  southern  M-dwarfs  only 
included  stars  with  V  <  14  mag  (Bonfils  et  al.  2013).  With  the 
ratio  of  early  M-dwarfs  (M0-M4)  to  mid  and  late  M-dwarfs 
(M4.5-M9)  being  23/28  for  stars  within  5  pc  (Cantrell  et  al. 
2013),  about  half  of  all  M-dwarfs  are  left  out  of  high-dispersion 
spectroscopic  surveys. 

To  address  the  question  of  M-dwarf  multiplicity  more  com¬ 
pletely,  we  initiated  a  volume-limited,  astrometric  and  infrared 
spectroscopic  search  for  companions  with  orbital  periods  be¬ 
tween  3  days  and  8  yr  around  mid  M-dwarfs  within  10  pc,  which 
includes  GJ  867B  as  a  target  (C.  L.  Davison  et  ah,  in  prepara¬ 
tion).  A  first  result  of  this  work  is  that  GJ  867B  is  revealed 
to  be  a  short-period  spectroscopic  binary  (SB).  We  discuss  the 
known  characteristics  of  the  GJ  867  system  in  Section  2.  We  list 
the  astrometric  results  in  Section  3  and  the  photometric  results 
in  Section  4,  and  newly  assigned  spectral  types  in  Section  5. 
We  discuss  the  spectroscopic  results  of  GJ  867B  and  the  new 
companion,  GJ  867D,  in  Section  6.  In  Section  7,  we  summarize 
the  results  and  discuss  the  case  for  physical  association  of  this 
hierarchical  system. 


6  LHS1070  was  claimed  to  be  a  quadruple  system  by  Henry  et  al.  (1999),  but 
later  determined  to  only  be  a  triple  system  (Kohler  et  al.  2012).  GJ  896A  and 
GJ  896B  were,  on  initial  inspection,  both  claimed  to  be  spectroscopic  binaries 
(Delfosse  et  al.  1999)  because  both  stars  showed  very  large  variations  in  radial 
velocity.  These  variations  were  later  thought  to  be  caused  by  inhomogeneous 
surface  features  combined  with  a  large  magnetic  spot  (X.  Delfosse  2013, 
private  communication). 


2.  THE  GJ  867  SYSTEM 

GJ  867A  (Vkc7 8  =  9.09  mag;  Bessel  1990)  is  die  primary 
component  of  a  widely  separated  visual  binary.  It  was  clas¬ 
sified  spectroscopically  by  Vyssotsky  &  Mateer  (1952)  as  an 
M-type  star;  we  more  accurately  classify  this  star  as  M2.0V 
(see  Section  5).  This  star  is  a  member  of  young  disk  (Ki- 
raga  &  Stepien  2007).  Its  visual  companion,  GJ  867B  ( Vkc  — 
11.45  mag;  Bessel  1990),  is  separated  by  24"5  at  a  position 
angle  (RA.)  of  350.45,  based  on  Two  Micron  All  Sky  Survey 
position  measurements.  The  spectral  type  of  GJ  867B  is  M3.5V 
(see  Section  5).  Coordinates,  magnitudes,  and  stellar  properties 
of  these  stars  are  given  in  Table  2. 

The  two  stars  are  presumed  to  be  physically  associated  from 
their  large  (~0'.'45  yr-1)  and  very  similar  proper  motions  (Dyer 
1954;  van  Altena  et  al.  1995;  Roser  et  al.  2011;  Zacharias 
et  al.  2013).  Nevertheless,  the  stars  have  slightly  but  steadily 
changed  their  relative  position  since  first  recorded  in  1 877  with 
a  separation  of  21"4  and  a  P.A.  of  359°,  which  presumably 
represents  the  orbital  motion  of  the  system  (Herbig  &  Moorhead 
1965). 

GJ  867A  was  discovered  by  Herbig  &  Moorhead  (1965)  to  be 
a  double-lined  SB  with  an  orbital  period  of  4.083  days.  The  two 
components  (A  and  C)  have  an  estimated  mass  ratio  ( MC/MA ) 
of  0.80  based  on  their  relative  velocity  amplitudes.  GJ  867B 
was  identified  with  high-resolution  optical  spectroscopy  as  a 
candidate  SB  by  Gizis  et  al.  (2002),  and  later  was  listed  as  either 
having  a  high  rotational  velocity  or  being  an  SB  by  Bonfils  et  al. 
(2013).  However,  until  now,  its  status  as  a  short-period  binary 
remained  unconfirmed. 

3.  ASTROMETRY 
3.1.  Observations 

Direct  images  of  GJ  867BDS  were  obtained  at  the  0.9  m 
telescope  at  CTIO  from  2003  July  until  2012  September, 
yielding  a  total  of  16  astrometric  epochs  over  9.17  yr.  Images 
were  obtained  using  the  2048  x  2048  Tektronix  CCD  camera 
on  the  central  quarter  of  the  CCD  chip.  The  CCD  chip  has 
a  plate  scale  of  O'.' 401  pixel  ,  which  gives  a  field  of  view 
of  6(8  x  6' 8  (Jao  et  al.  2003).  Astrometric  measurements  of 
GJ  867BD  were  made  in  the  V j  filter,  which  provided  a  field 
of  five  suitable  reference  stars.  No  astrometric  measurements 
were  taken  of  GJ  867 AC,  as  this  star  saturates  before  sufficient 
signal-to-noise  ratios  (S/Ns)  are  obtained  for  the  reference  stars . 
Each  night  that  GJ  867BD  was  observed,  we  obtained  5  to 
10  images,  for  a  total  of  96  images  on  16  different  nights. 
Exposure  times  for  an  individual  image  varied  from  25  to  100  s, 
with  an  average  time  of  54  s.  We  aimed  to  observe  all  images 
within  30  minutes  of  the  star  crossing  the  meridian  to  minimize 
the  effects  of  differential  color  refraction  (Jao  et  al.  2005).  In 
each  image,  the  FWHM  of  GJ  867BD  is  less  than  2"4;  this 
size  is  typically  set  by  the  seeing.  The  ellipticity  of  the  target 
star  and  all  reference  stars  are  less  than  20%.  These  criteria  are 
adopted  to  eliminate  elongated  reference  stars  (e.g.,  binaries), 
and  to  avoid  data  with  tracking/guiding  errors. 

During  the  time  of  observations,  two  different  Vj  filters 
were  used.  The  first  V j  filter  possessed  a  small  crack  and  was 
replaced  in  2005  February  with  a  nearly  identical  Vj  filter. 


7  Subscripts:  KC  indicates  Kron-Cousins  (SAAO  system),  and  used  later  in 
the  text,  J  indicates  Johnson. 

8  When  the  combined  light  of  the  individual  components  of  the  binary  pair 
cannot  be  separated,  we  give  the  object  the  name  of  the  binary  pair. 
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Table  2 

Characteristics  of  the  GJ  867  System 


Star 

R.A. 

Decl. 

Vkc 

Rkc 

!  KC 

J 

H 

Ks 

Spectral 

v  sin  i 

(J2000.0) 

(J2000.0) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

Class 

(km  s_1) 

GJ  867 AC 

22  38  45.601 

-20  37  16.1 1 

9.092 

8.082 

6.882 

5.671 

5. 1 1 1 

4.801 

M2.0VJ3 

4.74 

GJ  867BD 

22  38  45.3 1 1 

-20  36  51.91 

1 1.476 

10.296 

8.776 

7.341 

6.821 

6.49* 

M3.5VJ3 

8.73 

References.  (1)  Cutri  et  al.  2003;  (2)  Bessel  1990;  (3)  This  work;  (4)  Houdebine  2010. 


which  we  confidently  adopted  knowing  that  it  possessed  the 
same  transmission  spectrum  as  the  original.  After  a  few  years 
of  obtaining  data,  we  noticed  a  few  milliarcsecond  (mas)  offset 
in  the  astrometric  residuals  of  some  stars  that  were  known  to 
be  single  from  other  techniques  (Subasavage  et  al.  2009);  we 
determined  that  the  offset  first  appeared  when  we  replaced  the 
V j  filter.  The  effect  of  using  the  replacement  V j  filter  is  more 
clearly  seen  in  the  astrometric  residuals  in  the  right  ascension 
axis.  After  a  close  inspection  of  the  first  V j  filter,  it  was  realized 
that  the  small  crack  does  not  project  onto  the  part  of  the  chip  that 
is  regularly  used  (i.e.,  the  central  quarter).  Therefore,  in  2009 
July,  we  replaced  the  second  V j  filter  with  the  original  V j  filter 
(Riedel  et  al.  2010).  All  data  acquired  after  2009  July  were  with 
the  original  V j  filter.  Using  data  from  both  filters  may  increase 
the  parallax  error,  but  does  not  change  the  parallax  value.  We 
use  all  points  in  this  analysis  to  maximize  the  parallactic  ellipse 
coverage. 

3.2.  Astrometric  Reductions 

Here,  we  briefly  summarize  the  reduction  technique;  further 
details  can  be  found  in  Jao  et  al.  (2005)  and  Subasavage 
et  al.  (2009).  We  obtain  photometric  parallaxes  using  the  CCD 
distance  relations  of  Henry  et  al.  (2004)  from  which  distances 
are  calculated  for  all  of  our  reference  stars.  We  require  reference 
stars  to  be  further  than  100  pc  from  the  Sun,  and  initially  assume 
the  sums  of  all  reference  stars’  proper  and  parallactic  motions 
to  be  zero  mas  yr-1.  We  measure  the  centroid  of  our  target  star 
against  the  presumed  fixed  grid  of  five  reference  stars  using  the 
SExtractor  centroiding  algorithm  described  in  Bertin  &  Arnouts 
(1996).  Relative  to  this  reference  grid,  we  obtain  the  proper 
motion,  parallax  and  any  residual  astrometric  motion  of  the 
target  star  using  the  Gaussfit  program  from  Jefferys  et  al.  (1987). 
We  use  the  photometric  parallaxes  of  the  reference  stars  from  our 
VRI  photometry  and  correct  for  the  parallactic  motions  of  our 
reference  stars  to  determine  the  absolute  parallax  of  GJ  867BD. 

3.3.  Parallax  and  Proper  Motion  of  GJ  867BD 

We  measure  a  relative  parallax  for  GJ  867BD  of  109.45  ± 
1.74  mas.  The  correction  for  parallactic  motion  is  0.93  ±  0.20, 
yielding  a  final  absolute  parallax  of  1 10.38  ±  1 .75  mas.  This  is 
the  first  parallax  measurement  for  this  star  and  corresponds 
to  a  trigonometric  distance  of  9.06  ±  0.15  pc.  The  errors 
on  the  parallax  measurement  represent  our  internal  statistical 
precision.  To  estimate  our  systematic  errors,  we  determine  the 
standard  deviation  of  the  differences  between  the  parallaxes  of 
previous  RECONS  measurements  (e.g.,  Subasavage  et  al.  2009; 
Riedel  et  al.  2010;  Jao  et  al.  2011)  and  HIPPARCOS,  which  is 
4.3  mas.  This  error  dominates  our  internal  error  and  we  use  it 
here  to  compare  our  values  to  previously  published  values.  Two 
parallax  values  are  available  for  GJ  867 AC,  1 15. 10  ±  7.40  mas 
(van  Altena  et  al.  1995)  and  115.01  ±  1.32  mas  by  HIPPARCOS 
(van  Leeuwen  2007).  Our  parallax  measurement  of  GJ  867BD 
is  within  1.1  o  of  the  weighted  parallax  measurement  from  van 


Altena  et  al.  (1995)  and  van  Leeuwen  (2007)  for  GJ  867 AC. 
The  similar  distance  measurements  strengthen  the  case  that 
GJ  867AC  and  GJ  867BD  are  physically  associated.  The 
weighted  mean  parallax  of  the  GJ  867  system  is  113.37  ± 
1 .04  mas,  which  corresponds  to  a  distance  of  8.82  ±  0.08  pc. 

Our  measurements  of  the  proper  motion  amplitude,  p,  and 
the  proper  motion  RA.  for  GJ  867BD  are  421.8  ±  0.6  mas  yr~' 
and  97?  1  ±  0?1,  respectively.  The  errors  on  the  proper  motion 
amplitude  represent  our  internal  statistical  precision.  To  estimate 
our  systematic  errors,  we  again  use  the  standard  deviation 
of  the  differences  between  the  proper  motion  amplitudes  of 
previous  RECONS  measurements  (e.g.,  Subasavage  et  al.  2009; 
Riedel  et  al.  2010;  Jao  et  al.  2011)  and  HIPPARCOS ,  which  is 
16.8  mas  yr-1.  Similar  to  our  parallax  errors,  the  systematic 
error  for  the  proper  motion  amplitude  dominates  the  error 
budget.  Salim  &  Gould  (2003)  measure  the  ji  of  GJ  867BD 
to  be  431.2  ±  10.2  mas  yr-1  and  the  proper  motion  P.A.  to 
be  98?  1,  which  is  within  0.6a  of  our  measurement  of  the  p  of 
GJ  867BD;  they  do  not  provide  an  uncertainty  on  their  measured 
position  angle.  GJ  867AC  has  a  p  of  455.8  ±  2.7  mas  yr-1 
and  a  position  angle  of  99?9  (Salim  &  Gould  2003).  The 
difference  in  proper  motion  amplitudes  for  the  AC  and  BD 
components  is  ~35  mas  yr-1,  which  is  slightly  larger  than  our 
systematic  error  of  16.8  mas  yr-1  and  nevertheless  consistent 
with  orbital  motion.  Table  3  summarizes  previous  parallax  and 
proper  motion  measurements. 

3.4.  Detection  Limits  for  Astrometric  Companions 

The  precise  position  measurements  of  GJ  867BD  spanning 
almost  10  yr  allow  us  to  set  mass  limits  on  long-period  astro- 
metric  companions.  This  first  requires  an  accurate  assessment  of 
the  position  uncertainties  which  we  estimate  empirically.  This  is 
done  by  characterizing  the  statistical  uncertainty  associated  with 
each  epoch’s  position  measurement.  To  do  this,  we  calculate  the 
mean  position  error  per  night,  which  is  the  average  standard  de¬ 
viation  of  position  offsets,  once  parallactic  and  proper  motions 
are  removed.  The  mean  position  errors  for  all  nights  in  R.A.  and 
decl.  are  2.31  and  4.00  mas,  respectively.  The  average  resid¬ 
ual  deviation  from  the  position  predicted  by  the  proper  motion 
and  parallax  fit  is  calculated  as  the  standard  deviation  of  the 
nightly  offsets,  and  in  R.A.  and  decl.  are  1.37  and  3.83  mas, 
respectively.  As  we  are  measuring  the  predicted  position  with 
better  precision  than  the  average  error,  this  may  be  an  indicator 
that  we  are  slightly  overestimating  the  mean  position  errors.  See 
Figure  1. 

In  order  to  set  limits  on  the  presence  of  companions  orbiting 
GJ  867BD  with  orbital  periods  between  2  to  10  yr,  we  perform 
100,000  Monte  Carlo  simulations  with  a  range  of  orbital 
properties.  Orbital  periods  are  restricted  to  2,  4,  6,  8,  and 
10  yr,  which  are  representative  of  our  observational  sampling 
and  baseline.  We  estimate  the  mass  of  the  primary  to  be 
0.29  ±  0.06  Mq  based  on  the  mass  luminosity  relation  of  Henry 
&  McCarthy  (1993).  Although  we  adopt  this  mass  for  all  of  our 
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Table  3 

Parallax  and  Space  Motions  of  the  GJ  867  System 


Group 

Targets 

Abs.  7T 
(mas) 

(mas  yr-1) 

PA. 

(deg) 

Vean 

(kmC1) 

Reference 

RECONS 

GJ  867 AC 

RECONS 

GJ  867BD 

110.38  ±  1.75 

421.8  ±0.6 

97.1  ±0.1 

18.1 

This  paper 

HIP 

GJ  867 AC 

115.01  ±  1.32 

455.9  ±  1.7 

100.0 

18.8 

van  Leeuwen  (2007) 

HIP 

GJ  867BD 

NLTT 

GJ  867 AC 

455.8  ±  2.7 

99.9 

Salim  &  Gould  (2003) 

NLTT 

GJ  867BD 

431.2  ±  10.2 

98.1 

Salim  &  Gould  (2003) 

YPC 

GJ  867 AC 

115.10  ±7.40 

459 

97 

18.9 

van  Altena  et  al.  (1995) 

YPC 

GJ  867BD 

449 

99 

van  Altena  et  al.  (1995) 

LEP 

GJ  867 AC 

458.0 

100 

Lepine  &  Gaidos  (2011) 

LEP 

GJ  867BD 

Lepine  &  Gaidos  (201 1) 

TYCHO-2 

GJ  867 AC 

455.8  ±  2.7 

99.9 

H0g  et  al.  (2000) 

TYCHO-2 

GJ  867BD 

H0g  et  al.  (2000) 

UCAC4 

GJ  867 AC 

457.6  ±  11.3 

100.0 

Zacharias  et  al.  (2013) 

UCAC4 

GJ  867BD 

432.6  ±  15.2 

100.6 

Zacharias  et  al.  (2013) 

PPMXL 

GJ  867 AC 

454.9  ±  2.7 

100.2 

Roser  et  al.  (2011) 

PPMXL 

GJ  867BD 

419.1  ±4.6 

98.4 

Roser  et  al.  (201 1) 

Weighted  mean 

GJ  867 

113.37  ±  1.04 

426.2  ±  0.4 

99.2 

18.6 

Notes.  Several  numbers  related  to  characterizing  the  motions  of  this  system  are  given  without  errors.  All  parallax  and  proper  motion  values 
given  with  errors  of  GJ  867 AC  and  GJ  867BD  are  included  in  the  calculation  for  the  weighted  mean  values  of  GJ  867.  The  weighted  mean 
position  angle  and  weighted  mean  tangential  velocities  values  are  pure  averages  of  the  available  numbers,  as  none  of  these  numbers  have  listed 
errors. 
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Figure  1.  Astrometric  residuals  of  GJ  867BD  after  solving  for  the  proper  motion 
and  parallactic  motion.  The  mean  position  errors  in  R.A.  and  decl.  are  2.31  mas 
and  4.00  mas,  respectively.  These  are  slightly  larger  than  the  average  residual 
deviations  about  zero,  which  are  1.37  mas  and  3.83  mas  for  R.A.  and  decl., 
respectively.  From  these  data  we  see  no  indication  of  long  term  variations  and 
exclude  the  presence  of  most  brown  dwarf  and  massive  planets  with  long  orbital 
periods  (2-10  yr). 


calculations,  we  note  that  this  mass  is  based  on  the  brightness  of 
both  GJ  867B  and  GJ  867D,  as  we  are  unable  to  deconvolve  the 
binary  and  determine  magnitude  and  color  estimates  for  each 
component.  The  companion  star,  GJ  867D,  appears  to  contribute 
little  light  to  the  system  (see  Section  6.3),  and  thus  is  unlikely  to 
severely  bias  this  estimate.  This  mass  determination  is  slightly 
higher  than  the  mass  determined  by  the  spectral  type  of  GJ  867B 
and  the  spectral  type  mass  relations  in  Kraus  &  Hillenbrand 
(2007)  of  ~0.25  M0,  but  consistent  within  the  errors.  We  use  the 
weighted  distance  measurement  of  the  GJ  867  system  reported 
in  Section  3.3.  We  allow  the  inclination  to  vary  between  0°  and 
90°  and  the  eccentricity  to  vary  between  0  and  0.9.  We  allow  the 
longitude  of  periastron,  co,  and  the  longitude  of  the  ascending 
node,  £2,  to  vary  from  0°  to  180°.  This  simulation  suggests  that 


90%  of  the  time,  we  would  have  been  able  to  detect,  with  99.7% 
(or  3er )  confidence,  a  companion  with  a  mass  of  at  least  18  Mjup 
in  a  2  yr  period,  11.5  MJUP  in  a  4  yr  period,  9  MJUP  in  a  6  yr 
period,  7.5  Mjup  in  a  8  yr  period,  and  7  Mjup  in  a  10  yr  period. 
The  implication  is  that  GJ  867BD  is  unlikely  to  have  any  brown 
dwarf  or  massive  planetary  companions  with  periods  between  2 
to  10  yr.  The  orbital  periods  correspond  to  separations  of  1.1  to 
3. 1  AU  (0'.'07  to  0'.'22,  if  face  on). 

4.  NEW  SYSTEM  PHOTOMETRY 
4.1.  Photometric  Observations 

We  obtained  photometric  measurements  of  GJ  867BD  and 
its  reference  stars  on  three  different  nights  using  an  identical 
setup  as  that  used  for  the  astrometric  measurements  at  the 
CTIO  0.9  m  telescope.  As  our  initial  goal  was  to  obtain 
photometry  of  GJ  867BD,  we  allowed  GJ  867 AC  to  saturate  in 
our  images,  and  therefore  provide  no  new  apparent  magnitude 
measurements  for  GJ  867AC.  We  obtained  Vj.  Rkc,  and  7kc 
photometry  of  GJ  867BD  that  were  calibrated  using  at  least  10 
photometric  standard  stars  per  night  from  Landolt  (1992)  and 
Graham  (1982).  All  images  were  obtained  at  an  airmass  of  less 
than  2. 


4.2.  Photometric  Reductions 

Observations  of  the  photometric  standards  are  used  to  de¬ 
termine  transformation  equations  and  to  create  the  extinc¬ 
tion  curves,  enabling  us  to  calculate  absolute  photometry  for 
GJ  867BD.  Specific  details  are  given  in  Jao  et  al.  (2005)  and 
Winters  et  al.  (2011).  The  apparent  magnitudes  through  the  Vj, 
Rkc,  and  7kc  filters  are  11.47  mag,  10.29  mag,  and  8.77  mag, 
respectively,  with  an  error  of  0.01  mag  each.  Our  measurements 
are  within  0.02  mag  of  the  measured  photometry  by  Bessel 
(1990). 
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Although  astrometric  observations  do  not  allow  absolute 
photometry,  relative  photometry  can  be  extracted.  Consequently, 
all  images  used  to  search  for  variability  are  in  the  Vj  filter.  We 
measure  the  standard  deviation  of  the  instrumental  magnitude 
of  GJ  867BD,  (cr mag),  relative  to  reference  stars  for  over  16 
epochs  spanning  9.17  yr  to  assess  the  amount  of  photometric 
variability.  For  GJ  867BD,  we  use  the  same  set  of  reference  stars 
in  the  astrometric  reduction.  The  average  crmag  is  32.4  mmag, 
which  is  above  the  mean  activity  level  of  most  field  dwarfs 
(~13  mmag;  Jao  et  al.  2011),  but  is  not  surprising  given  its 
classification  as  a  flare  star  (Byrne  1978). 

5.  NEW  SPECTRAL  TYPES 

New  spectral  types  are  determined  using  the  RECONS 
classification  system,  which  focuses  on  K-  and  M-dwarfs.  The 
spectral  region  used  for  this  covers  6000-9500  A  at  a  resolution 
of  8.6  A;  further  details  are  given  in  Jao  et  al.  (2008)  and  A.  R. 
Riedel  et  al.  (in  preparation).  Observations  were  made  using  the 
R-C  Spectrograph  on  the  CTIO  1.5  m  with  the  Loral  1200  x 
800  CCD  camera,  utilizing  the  #32  grating  in  first  order  with 
tilt  15°  1,  and  order  blocking  filter  OG570.  We  remove  the  Ha 
feature  and  the  telluric  features  on  the  spectra  using  the  sky 
transmission  map  from  Hinkle  et  al.  (2003).  We  then  compare 
our  stellar  spectra  to  our  library  of  standards  spectra  from  Henry 
et  al.  (2002)  and  determine  that  the  best  match  with  the  lowest 
standard  deviation  to  be  the  spectral  type  of  our  star.  Spectral 
types  are  accurate  to  0.5  spectral  subclasses.  We  determine 
the  spectral  types  of  GJ  867 AC  and  GJ  867BD  based  on  the 
combined  light  of  the  spatially  unresolved  pairs  to  be  M2.0  and 
M3. 5,  respectively.  These  values  are  consistent  to  within  0.5 
subclasses  of  previously  adopted  values  (e.g.,  Reid  et  al.  2004). 

6.  HIGH-DISPERSION  SPECTROSCOPY 

6.1.  Observations 

Spectrographic  observations  were  made  at  NASA’s  Infrared 
Telescope  Facility,  a  3  m  telescope,  using  the  high-resolution 
infrared  echelle  spectrograph,  CSHELL  (Tokunaga  et  al.  1990; 
Greene  et  al.  1993).  Spectra  were  centered  at  2.298  /im 
(vacuum)  and  span  approximately  50  A.  This  wavelength  region 
was  chosen  to  obtain  precise  radial  velocities  (RVs)  by  measur¬ 
ing  the  positions  of  the  rich  photospheric  12CO  R  lines  (2.3  //  m ) 
relative  to  a  wavelength  solution  set  by  telluric  methane  absorp¬ 
tion  features  in  the  Earth’s  atmosphere  (e.g.,  Blake  et  al.  2010; 
Bailey  et  al.  2012).  Observations  were  made  using  the  0/15  slit, 
yielding  a  measured  resolving  power  of  ~46,000.  Ten  measure¬ 
ments  of  GJ  867B  were  obtained  spanning  a  temporal  baseline 
of  100  days  between  2012  August  22  and  2012  November  29. 
Each  measurement  consists  of  two  spectra  taken  in  quick  suc¬ 
cession  of  one  another  along  the  vertical  direction  of  the  slit, 
separated  by  10".  The  two  different  positions  are  referred  to  as 
nod  positions.  The  first  three  RV  measurements  were  obtained 
during  late  2012  August,  over  a  4  day  period.  At  that  time,  we 
noticed  a  large  RV  variation  of  almost  18  km  s-1  and  flagged 
this  variable  candidate  as  a  high  priority  target  for  our  upcoming 
runs  to  confirm  it  as  an  SB.  Exposure  times  averaging  ~500  s 
each  were  set  to  obtain  predicted  S/Ns  of  ~100  per  exposure, 
and  resulted  in  pair-combined  S/Ns  ranging  from  70  to  230. 
Along  with  observations  of  GJ  867B,  spectra  of  either  the  A7V 
star  Altair  or  the  Al  V  star  Sirius  were  obtained  to  gain  an  initial 
estimate  of  the  wavelength  solution  and  the  instrumental  profile 
of  CSHELL,  since  these  spectra  contain  only  telluric  features. 


Finally,  20  flats  and  darks  with  integration  times  of  10  or  20  s 
were  taken  each  night. 

6.2.  Spectroscopic  Reductions 

A  brief  overview  highlighting  the  aspects  of  our  spectral 
extraction  and  RV  fitting  routine  will  be  given  here.  For  a  full 
explanation  of  the  data  reduction  technique,  the  reader  is  referred 
to  Bailey  et  al.  (2012)  and  C.  L.  Davison  et  al.  (in  preparation). 

To  prepare  the  two-dimensional  images  for  spectral  extrac¬ 
tion,  we  first  correct  for  flat  fielding  using  dark-subtracted, 
median-combined  flat  field  images.  Then,  we  remove  sky  emis¬ 
sion,  detector  bias,  and  dark  current  from  the  spectral  images  by 
subtracting  each  individual  image  from  its  pair  image  taken  at 
a  different  nod  position  that  immediately  preceded  or  followed 
the  image.  We  assume  that  the  changes  in  the  detector  and  spec¬ 
trograph  properties  are  negligible  over  the  short  period  in  which 
the  two  spectra  were  taken  at  the  different  nod  positions.  Then, 
we  perform  an  optimal  extraction  of  the  spectra  as  implemented 
for  nod  subtracted  pairs  (Horne  1986;  Bailey  et  al.  2012;  Tanner 
et  al.  2012).  Unlike  a  standard  extraction  that  simply  sums  over 
all  the  pixels  to  determine  the  one-dimensional  spectral  profile, 
an  optimal  extraction  sums  pixels  weighted  by  a  variance  image 
of  the  spectrum’s  spatial  profile  to  determine  the  spectral  pro¬ 
file.  It  is  assumed  that  the  spectral  profile  varies  smoothly  in  the 
direction  parallel  to  the  dispersion.  Using  a  weighted  image  to 
determine  the  spectral  profile  gives  the  advantage  that  it  mini¬ 
mizes  the  noisy  contribution  of  the  profile  wings  and  eliminates 
the  majority  of  the  features  on  the  spectral  profile  caused  by 
cosmic  rays  or  hot  pixels. 

6.3.  Spectroscopic  Analysis  and  Modeling  as  a 
Single-lined  Spectroscopic  Binary 

If  GJ  867BD  is  a  double-lined  SB,  we  would  expect  it  to 
appear  double-lined  in  at  least  some  of  our  spectra.  Visual 
inspection  showed  no  sign  of  a  secondary  component.  Likewise, 
if  its  companion  is  of  comparable  brightness,  we  might  observe 
broader  absorption  lines  when  the  stars  in  the  system  are  at 
opposition.  Broader  absorption  lines  would  then  be  translated 
by  our  analysis  code  into  larger  v  sin  i  values.  However,  all 
10  v  sin  i  measurements  are  consistent  to  within  3.1  km  s'1 
and  show  no  apparent  correlated  variation  with  the  changes 
in  RV.  These  factors  lead  us  to  believe  that  this  star  is  a 
single-lined  spectroscopic  binary  (SB1).  Also,  there  is  a  lack  of 
obvious  remaining  residuals  in  the  single  star  spectral  modeling, 
described  below. 

We  subsequently  analyze  the  system  as  an  SB1.  After  ex¬ 
tracting  the  spectra,  we  fit  each  observation  to  high-resolution 
spectral  models  that  are  convolved  to  the  resolution  of  CSHELL. 
The  models  are  formed  by  combining  synthetic  spectra  created 
from  NextGen  (Hauschildt  et  al.  1999)  with  empirical  telluric 
models  extracted  from  ultra  high-resolution  KPNO/FTS  telluric 
spectra  of  the  Sun  (Livingston  &  Wallace  1991).  Our  synthetic 
model  is  composed  of  spectral  features  from  a  single  M3.5V 
star.  In  order  to  make  a  realistic  model,  we  assume  the  starting 
stellar  parameters  of  GJ  867B  as  follows:  We  estimate  the  stellar 
temperature  to  be  3200  K  based  on  its  spectral  type  of  M3.5V, 
using  the  relationship  from  Kraus  &  Hillenbrand  (2007);  we 
set  its  surface  gravity  (log(g))  to  be  4.8  (cgs),  as  this  number  is 
consistent  with  other  field  dwarfs  of  similar  mass  (Hillenbrand 
&  White  2004).  Our  model  spectrum  consists  of  19  parameters. 
The  limb  darkening  coefficient  parameter  is  set  to  0.6,  a  reason¬ 
able  assumption  for  cool  stars  at  infrared  wavelengths  (Claret 
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Figure  2.  Spectra  of  GJ  867BD.  Spectra  are  modeled  by  combining  a  telluric 
spectrum  (top  spectrum),  which  provides  an  absolute  wavelength  reference, 
with  a  theoretical  stellar  spectrum  (second  spectrum).  The  CSHELL  spectrum 
of  GJ  867B  is  shown  (black;  third  spectrum)  in  comparison  with  the  best  fit  (red; 
third  spectrum).  The  mean  residuals  of  the  best-fit  spectrum  are  1.45%  (bottom 
spectrum).  The  average  mean  residuals  for  all  CSHELL  spectra  obtained  is 
3.09%. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

2000)  and  the  remaining  1 8  parameters  are  left  to  vary.  Of  those 
variable  parameters,  3  of  those  parameters  are  used  to  solve  the 
wavelength  solution  and  9  of  those  parameters  are  Gaussians 
used  to  model  the  point-spread  function  of  the  spectrum.  The 
remaining  6  are  the  depth  of  the  telluric  spectra,  depth  of  the 
synthetic  spectra,  v  sin  i,  RV,  and  2  continuum  normalization 
constants. 

To  determine  the  12  parameters  related  to  the  wavelength  and 
instrumental  profile,  we  first  fit  the  empirical  telluric  spectrum 
to  A  star  observations.  We  obtain  the  optimal  parameters  by 
minimizing  the  variance  weighted  reduced  /  2.  On  the  one  night 
when  an  A  star  was  not  observed,  we  averaged  parameters 
determined  on  nights  with  A  star  observations  and  during  the 
same  run  together  to  estimate  initial  guesses  for  the  wavelength 
solution  and  instrumental  profile.  This  step  is  critical,  as  our 
optimization  routine  is  sensitive  to  initial  guesses. 

After  obtaining  these  initial  parameters,  we  fit  the  combined 
model  of  the  telluric  and  synthetic  spectra  to  our  observed 
CSHELL  spectra.  During  this  iterative  process,  we  first  fit  for 
the  wavelength  solution,  the  depth  of  the  telluric  and  synthetic 
spectra,  the  RV,  and  the  two  normalization  constants.  Then,  on 
our  second  iteration,  we  fit  for  v  sin ;,  the  depth  of  the  telluric  and 
model  spectral  features,  and  the  two  normalization  constants. 
On  our  final  iteration,  we  optimize  the  wavelength  solution, 
the  depth  of  the  telluric  and  synthetic  spectra,  the  RV,  and  the 
two  normalization  constants.  All  of  our  fits  are  optimized  using 
AMOEBA,  which  uses  a  downhill  simplex  method  to  find  the 
best  solution  for  multiple  variables  (Nelder  &  Mead  1965).  We 
restrict  AMOEBA  to  specified  ranges  of  physically  reasonable 
solutions.  An  example  of  an  optimally  extracted  spectrum  of 
GJ  867BD,  the  telluric  model,  the  stellar  model,  and  the  best-fit 
combination  are  shown  in  Figure  2. 

6.4.  Spectroscopic  Results 

The  determined  RVs  and  projected  rotational  velocities 
( v  sin  i  values)  are  listed  in  Table  4.  The  observed  RV  error 
is  a  combination  of  theoretical  photon  noise  error,  intrinsic  stel¬ 
lar  error,  and  instrumental  error  (er0bs2  =  cr photon2  +  cr  stellar2  + 
cr instr2)-  The  photon  error  is  calculated  based  on  the  prescription 
by  Butler  et  al.  (1996).  We  assume  the  intrinsic  stellar  error  to 
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Figure  3.  RV  measurements  of  GJ  867B  vs.  barycentric  Julian  date.  A  Keplerian 
curve  is  fit  to  the  RV  measurements,  centered  on  the  systemic  velocity  of 
—6.1  km  s-1.  The  errors  on  the  RV  measurements  are  smaller  than  the  plot 
symbols. 


Table  4 

Spectroscopic  Results  for  GJ  867BD 


JD-240000 

Orbital 

Phase 

Radial  Velocity 
(ms-1) 

S/N 

v  sin  i 
(km  s_1) 

56161.36 

1.00 

7594  ±  49 

177 

9.3 

56163.41 

1.25 

-10132  ±53 

168 

8.4 

56164.45 

0.49 

7645  ±  61 

135 

9.8 

56173.32 

0.40 

1396  ±  53 

162 

7.6 

56173.41 

0.49 

7274  ±  62 

138 

8.9 

56174.31 

1.39 

-19888  ±46 

204 

9.3 

56174.41 

1.48 

-24330  ±  44 

228 

8.4 

56259.26 

0.22 

-14613  ±47 

127 

9.5 

56260.23 

1.20 

-3936  ±41 

197 

9.5 

56261.32 

0.50 

6415  ±  139 

71 

6.7 

Note.  The  reported  RV  measurements  have  been  corrected  for  the  barycentric 
motion  of  Earth. 


be  zero.  We  know  this  is  an  underestimate  of  the  stellar  error  for 
GJ  867B,  which  is  an  active  flare  star.  However,  we  do  not  have 
enough  data  to  accurately  assess  the  error  caused  by  intrinsic 
stellar  error  for  this  star.  We  adopt  the  average  instrumental  error 
of  68  m  s_1  determined  for  10  stars  with  low  rotational  veloc¬ 
ities  from  our  survey  of  mid  M-dwarfs  (C.  L.  Davison  et  al., 
in  preparation)  as  the  instrumental  error  for  our  RV  measure¬ 
ments  of  GJ  867BD.  The  average  v  sin  i  value  from  our  fitting 
prescription  is  8.74  ±  0.98  km  s-1.  The  error  on  the  v  sin  i  mea¬ 
surement  is  set  from  the  standard  deviation  of  the  independent 
measurements  from  all  10  epochs. 

6.5.  Periodicity  Analysis 
6.5.1.  The  Orbital  Period 

GJ  867B  shows  a  large  change  in  RV  on  timescales  of  less 
than  1  day  (Figure  3).  To  determine  if  the  changes  in  RV  are 
periodic,  we  use  a  Lomb-Scargle  algorithm  to  search  for  periods 
between  0.5  and  10  days.  The  highest  power  of  the  periodogram 
favors  a  period  of  1.795  days.  Unfortunately,  because  of  the 
limited  sampling,  there  are  many  other  peaks  of  comparable 
power  (Figure  4).  The  highest  peak  at  1.795  days  has  a  power 
value  of  4.37,  which  is  closely  followed  by  power  values  of 
4.24  and  4.21  for  the  next  most  likely  peaks  at  periods  of  1.541 
and  2.841  days,  respectively.  To  confirm  that  the  highest  peak  is 
representative  of  the  best-fit  period,  we  conduct  an  independent 
X1  analysis,  where  we  fit  circular  orbits  to  the  top  three  periods 
given  from  the  periodogram  fit.  Knowing  the  orbital  period 
is  less  than  3  days,  we  assume  the  orbits  are  circularized 
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Table  5 

Characteristics  of  the  Spectroscopic  Binaries 


Davison  et  al. 
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Figure  4.  Periodograms  of  the  RV  measurements  of  GJ  867BD.  The  left 
panel  illustrates  the  periodogram  for  the  full  range  of  periods  considered, 
while  the  right  panel  is  the  portion  containing  the  most  significant  period  of 
1.795  ±0.017  days. 


Phase  Phase  Phase 

Figure  5.  Phased  RV  curves  of  GJ  867BD.  Shown  are  the  RV  curves  of 
GJ  867BD  phased  to  1.541,  1.795,  and  2.841  days,  which  are  the  three  periods 
with  the  highest  significance  from  the  periodogram  analysis.  The  black  dotted 
line  is  the  best-fit  Keplerian  model;  the  best-fit  residuals  are  shown  in  the  bottom 
panels.  The  period  of  1.795  days  is  strongly  favored  because  of  its  much  smaller 
residuals. 

due  to  tidal  effects  (Mayor  &  Mermilliod  1984).  We  optimize 
the  circular  fitted  curves  to  give  the  lowest  observed-minus- 
calculated  (O  —  C)  RV  value  using  the  AMOEBA  minimization 
routine  (Section  3.4);  the  results  are  shown  in  Figure  5.  The 
average  (<9— C)  value  for  1.795  days  is  0.58  km  s  1  and  is 
almost  five  times  smaller  than  the  (O  —  C)  values  for  the  periods 
of  2.841  days  (2.88  km  s-1)  and  1.541  days  (3.00  km  s-1). 
Similarly,  the  /2  value  for  1.795  days  is  27  and  38  times  smaller 
than  the  /2  values  at  1.54  days  and  2.841  days,  respectively. 
From  these  comparisons,  we  conclude  that  1.795  days  is  the 
orbital  period  of  the  BD  system. 

6.3.2.  The  Uncertainties  in  the  Orbit 

To  determine  the  orbital  properties  and  the  associated  un¬ 
certainties  of  the  orbit,  we  use  the  nonlinear  least  squares 
fitting  IDL  routine,  mpfitfun.pro  (Markwardt  2009)  to  model 
our  data  as  a  sinusoidal  curve.  The  amplitude  of  the  fit¬ 
ted  curve  is  21.4  ±  0.5  km  s-1,  the  velocity  of  the  system 
is  —6.1  ±  0.3  km  s-1,  and  the  epoch  of  conjunction  is  JD 
2,456,160.5611  ±0.0268. 

The  peak  with  the  most  likely  period  is  in  a  cluster  of  aliased 
peaks.  To  test  whether  the  neighboring  peaks  represent  rea¬ 
sonable  orbital  solutions,  we  again  fit  circular  orbits  to  the 
top  three  aliasing  peaks  in  closest  proximity  to  the  favorable 
period.  The  (O  —  C)  values  for  these  peaks  range  from 
3.16  km  s-1  to  5.56  km  s“  1  making  these  values  significantly 
larger  than  the  ( O  —  C )  values  from  the  top  three  likely  pe¬ 
riods  given  by  the  periodogram  fit.  We  conclude  that  neither 


GJ  867 ACa 

GJ  867BD 

Period  (days) 

4.08322  ±  0.00004 

1.795  ±  0.017 

y  (km  s_1) 

-8.7 

-6.1  ±0.3 

Ki  (km  s~*) 

46.8  ±  0.45 

21.4  ±0.5 

Ki  (km  s~*) 

58.1  ±0.46 

T0  (days) 

JD  2437144  ±  0.006 

JD  2456160.5611  ±0.0268 

Mi  (M0) 

0.271  ±0.006/ sin3; 

0.29  ±  0.59b 

M2  (Mq) 

0.218  ±().005/sin3r 

>  0.056  ±  0.007b 

01  sin/  (km) 

2.63  (±0.03)  x  106 

0.528  (±  0.02)  x  106 

02  sin/  (km) 

3.26  (±0.03)  x  106 

e 

0.010  ±0.010 

0C 

O) 

356°  ±  32° 

0C 

Notes. 

a  Binary  properties  for  GJ  867 AC  are  from  Herbig  &  Moorhead  (1965). 
b  The  mass  estimate  of  the  primary  is  based  on  a  mass  luminosity  relation  while 
that  of  the  secondary  is  a  minimum  mass  (m  sin  i;  see  the  text). 
c  The  eccentricity  and  longitude  of  periastron  are  assumed  to  be  zero. 

of  these  periods  represent  good  fits  to  our  data  and  determine 
the  uncertainty  in  the  period  from  the  FWHM  of  the  central 
peak  only,  which  yields  a  final  best-fit  period  and  uncertainty  of 
1.795  ±0.017  days. 

The  average  difference  between  the  {()  —  C )  value  for  our 
best-fit  circular  orbit  of  0.58  km  s~  1  is  significantly  higher  than 
our  average  estimated  error  for  this  star  based  on  the  theoretical 
photon  noise  error  and  the  instrumental  error.  This  could  be 
caused  by  a  few  possibilities.  GJ  867B  could  have  a  slightly 
elliptical  orbit,  and  therefore  the  deviation  we  see  from  the 
observed  compared  to  calculated  value  is  real.  Alternatively, 
our  estimated  errors  could  be  low,  as  we  did  not  account  for  any 
intrinsic  stellar  error.  GJ  867B  is  a  flare  star  with  a  nominally 
higher  rotational  velocity,  which  could  be  an  indicator  of  a 
higher  amount  of  intrinsic  stellar  jitter  compared  to  other  older 
field  dwarfs.  We  caution  that  if  the  RV  errors  are  too  low,  then  the 
calculated  errors  on  the  orbital  properties  will  be  underestimated 
as  well. 

6.6.  Limits  on  the  Mass  of  GJ  867D 

Without  a  direct  detection  of  GJ  867D,  we  can  only  provide  a 
lower  mass  limit  based  on  the  observed  reflex  motion  of  GJ  867B 
using  Kepler’s  third  law.  We  assume  the  orbit  is  circular  and  that 
the  mass  of  the  primary  star  is  0.29  ±  0.06  M0  (Section  3.4). 
Under  these  assumptions  the  minimum  mass  of  GJ  867D  is 
61  ±  7  Mjup  (0.059  ±  0.007  M0);  the  uncertainty  in  the 
minimum  mass  is  determined  primarily  by  the  uncertainty  in 
the  primary  star’s  mass.  GJ  867D  could  be  a  brown  dwarf.  The 
spectral  properties  of  GJ  867 AC  and  GJ  867BD  are  summarized 
in  Table  5. 

The  brightness  of  GJ  867BD  supports  the  conclusion  that 
GJ  867D  is  very  low  mass.  We  compare  the  absolute  magnitude 
in  V  to  the  (V—K)  color  of  stars  within  5  pc  (Cantrell  et  al. 
2013)  to  GJ  867  system  (see  Figure  6).  GJ  867BD  lies  within  the 
broad  band  making  up  the  main  sequence.  Therefore,  there  is  no 
evidence  that  the  companion  contributes  much  light,  consistent 
with  its  undetected  spectral  features. 

7.  SUMMARY 

We  report  infrared  spectroscopic  observations  obtained  with 
the  CSHELL  spectrograph  that  reveal  GJ  867B  to  be  an  SB1 
with  an  orbital  period  of  1.795 ±0.017  days;  the  components  are 
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Figure  6.  Absolute  V  magnitude  vs.  (V—K)  colors.  This  color-absolute 
magnitude  diagram  shows  GJ  867 AC  and  GJ  867BD  relative  to  all  known 
M-dwarfs  within  5  pc.  This  plot  shows  that  GJ  867BD  is  not  significantly  over- 
luminous  relative  to  this  empirical  main  sequence,  consistent  with  GJ  867BD 
being  a  single-lined  SB  and  GJ  867D  being  low  mass. 

labeled  B  and  D.  Assuming  that  GJ  867B  has  a  mass  of  0.29  Mq, 
then  GJ  867D  has  a  minimum  mass  (m  sin  i)  of  61  ±  7  Mjup; 
it  could  be  a  brown  dwarf.  Astrometric  measurements  spanning 
nearly  10  yr  are  used  to  determine  the  trigonometric  parallax 
and  proper  motion  of  GJ  867BD.  The  trigonometric  parallax  of 
GJ  867BD  is  consistent  to  within  1.1a  of  the  measured  parallax 
of  GJ  867 AC  by  HIPPARCOS,  confirming  that  the  components 
of  this  24" 5  pair  are  at  similar  distances.  The  measured  proper 
motion  amplitude  of  GJ  867BD  is  421.8  ±  0.6  mas  yr-1, 
which  is  similar  to  but  slightly  smaller  than  the  proper  motion 
amplitude  of  GJ  867 AC  (455.8  ±  2.7  mas  yr-1;  Salim  &  Gould 
2003);  both  stars  are  moving  in  the  same  direction  to  within 
2°.  The  difference  (-34.0  mas  yr-1)  can  easily  be  explained 
by  orbital  motion.  Adopting  the  minimum  masses  in  Table  5 
and  assuming  the  projected  separation  (216  AU)  is  the  semi¬ 
major  axis  of  the  AC  and  BD  pair,  these  components  would  be 
moving  relative  to  one  another  at  1.8  km  s_1,  corresponding 
to  43.2  mas  yr-1  (at  8.82  pc)  if  in  a  face-on  circular  orbit.  If 
in  a  somewhat  inclined  orbit  (and  if  the  components  are  more 
massive  than  the  minimum  masses  used  here),  orbital  motion 
could  also  explain  the  difference  in  the  systemic  velocities  of 
GJ  867BD,  at  -6.1  ±0.3  km  s”1  (Table  5),  and  GJ  867 AC, 
at  —8.7  km  s_I  (Herbig  &  Moorhead  1965);  we  consider 
this  apparent  discrepancy  less  relevant,  however,  since  the 
systemic  velocity  for  GJ  867AC  is  based  on  much  coarser 
precision  individual  measurements  (1.4  km  s-1)  and  no  formal 
uncertainty  on  the  systemic  velocity  is  provided.  Altogether, 
the  similar  proximity,  distance,  and  space  motion  of  these 
two  high  proper  motion  SBs  establish  them  as  very  likely 
physically  associated.  As  such,  the  four  star  system  GJ  867 
becomes  one  of  only  four  quadruple  systems  within  10  pc, 
and  the  closest  quadruple  system  to  the  Sun  with  an  M-dwarf 
primary. 

The  discovery  of  this  nearby  quadruple  system  suggests 
that  many  more  hierarchical  low-mass  systems  may  exist 
nearby,  but  have  not  been  found  because  of  observational 
biases.  As  highlighted  in  the  introduction,  the  discovery  and 
characterization  of  systems  like  this  potentially  offer  a  powerful 
tool  for  constraining  the  uncertain  physics  involved  in  binary 
and  multiple  star  formation.  We  note  that  nearby  systems  like 
GJ  867  are  especially  valuable,  since  these  stars  are  close  enough 
to  have  the  spectroscopic  orbits  that  can  be  spatially  resolved 
interferometrically,  allowing  accurate  masses  of  these  systems 
to  be  determined. 
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